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ABSTRACT 


A  computer  simulation  of  a  pulsed  HF  laser  pumped  by  the  +  F^ 
chain  reaction  is  presented.  A  chemical  kinetic  model  encompassing  68 
reactions  is  used  to  approximate  the  reacting  mixture  contained  within  an 
optical  cavity.  For  each  vibrational  level,  a  Boltzmann  distribution  for  the 
rotational  levels  is  assumed  with  lasing  on  the  vibrational  band  at  line  center 
of  the  transition  having  maximum  gain. 

An  analysis  of  cavity  and  chemical  mechanisms  yields  a  simple  relation* 
ship  between  the  pumping  and  depletion  rates  of  the  vibrational  levels.  This 
relationship  is  used  to  make  a  comprehensive  study  of  the  effects  of  cavity 
and  chemical  parameters  on  the  laser  pulse.  The  effect  of  changes  in  un¬ 
certain  chemical  reaction  rate  coefficients  is  assessed.  A  competition  exists 
between  vibration-translation  deactivation  reactions  and  stimulated  emission. 
Laser  performance  is  most  sensitive  to  the  relative  rates  of  the  vibrational- 
translational  and  pumping  reactions.  In  addition,  the  effect  of  the  photon  flux 
on  the  chemical  mechanisms  is  significant. 
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NOMENCLATURE 
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EL(t> 


Jv+l,v 


k  ,  k 
n  fr 

k  >  k, 
-n  br 


Lfr,Lbi 


Einstein  isotrupic  absorption  coefficient 

speed  of  light,  2,  9979  X  lO^cm/sec 

specific  heat  at  constant  pressure  of  species  i 

chemical  deactivation  rate  from  level  v  +  1  to  v 

energy  emitted  during  laser  pulse  to  time  t 

rotational  energy  associated  with  vibrational  level 
v  and  rotational  level  J 

total  energy  emitted  during  lasing  on  a  v  +  1  — ►  v  band 
Planck's  constant,  6.  6256  X  10  joule-sec 
molar  enthalpy  of  species  i 
rotational  quantum  number 

rotational  quantum  number  associated  with  the  lower  level 
of  a  lasing  transition 

Boltzmann's  constant,  1.3805  X  *  0*  joule /°K 
forward  rate  coefficients 
backward  rate  coefficients 
length  of  active  medium 

forward  and  backward  rates  for  reaction  r 
catalytic  species  of  type  j 
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N, 


N. 
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N(v) 

N(v,  J) 


Pi 


v+l,v 

Q^(T) 


R0’  RL 


tl/2 


T. 

1 


th 


a  .,P  . 
ri  ri 


u(v,  J) 
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Avogadro's  number,  6.  0225  X  10  molecules / mole 
molar  concentration  of  species  i 

molar  concentration  of  lasing  species  with  vibrational 
level  v 

molar  concentration  of  lasing  species  with  vibrational 
level  v  and  rotational  level  J 

pressure 

initial  pressure  of  the  reacting  mixture 
output  power 

output  power  for  lasing  into  level  v 

rotational  partition  function  for  vibrational  level  v 

mirror  reflectivities 

time 

time  to  emit  first  1/2  of  the  total  pulse  energy- 

pulse  duration 

temperature 

initial  temperature  of  the  reacting  mixture 

vibrational  quantum  number 

threshold  gain 

stoichiometric  coefficients 

optical  gain  for  transitions  to  state  (v,  J) 


-x  - 


“n 


K 

V 


0(v,  J) 


4-rad(v'J) 
<*>c(v.  J) 


laser  efficiency 

rate  of  pumping  into  vibrational  level  v 

reaction  rate  weighting  factor 

Voigt  profile  for  level  (v,J) 

net  chemical  production  rate  (mole/cm  -sec) 
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I.  INTRODUCTION 


Increased  interest  in  conversion  of  chemical  energy  into  laser  power 

and  the  scientific  value  of  quantum-chemistry  data  available  from  the  use  of 

1-4  4-7 

the  chemical  laser  have  stimulated  recent  experimental  and  theoretical 
studies.  A  current  review  of  chemical  laser  technology  has  been  made  by 

g 

Dzhidzhov,  Platonenko,  and  Khokhlov.  The  multiplicity  of  phenomena 
operating  in  a  chemical  laser  makes  quantitative  predictions  of  laser  perfor¬ 
mance  or  detailed  interpretation  of  experimental  data  difficult  without  the 
aid  of  comprehensive  modeling  of  these  phenomena.  Early  modeling  studies 
estimated  the  performance  of  a  pulsed  laser  by  considering  the  time  history 

of  gain  on  a  given  band,  while  the  effect  of  the  cavity  and  induced  emission 

6  7 

on  the  chemistry  were  neglected.  '  More  recently,  computer  simulations 

4  9 

have  included  these  effects,  ’  and  such  a  code  is  used  in  this  study.  The 
simulation  is  used  to  study  the  effects  of  various  parameters  on  the  behavior 
of  a  pulsed  HF  laser.  The  principal  objectives  of  this  study  are  to 
determine: 

1.  Relative  importance  of  the  chemical  kinetic  mechanisms 
involved  in  the  pulsed  HF  laser 

2.  Effect  of  the  optical  cavity  on  the  reacting  +  F2  mixture 

3.  Effect  of  various  parameters  such  as  the  initial  concentrations 
and  the  initial  temperature. 

4 

The  method  of  simulation  is  similar  to  that  used  by  Airey  for  an  HCI 
laser.  The  major  provisions  of  the  laser  model  are  summarized  as  follows: 

i.  The  dominant  reactions  used  to  represent  the  chemical  processes 
in  the  reacting  mixture  are: 

a.  the  H2  -  F2  chain 

F  +  H2  ~ — (vl  +  H 

H  +  F2  ?=SHF(v)  +  F 


b.  vibrational-translational  (VT)  deactivation 

HF(v)  +  M  HF(v  -  1)  +  M 

c.  vibrational- vibrational  (VV)  quantum  exchange 

HF(v)  +  HF(v'  )  HF(v  +  1  )  +  HF(v'  -  1  ) 

HF (v)  +  H2(v'  )  ^==£  HF(v  +  1)  +  H2(v'  -  1) 

d.  dissociation-recombination 

F2  +  M  M  +  F  +  F 

H2  +  M  - - -  M  +  H  t  H 

HF(v)  +  M  ^:~=-  M  +  H  +  F 

This  reaction  system  is  approximated  by  the  68  reactions  listed 
in  the  Appendix. 

2.  The  reacting  mixture  is  homogeneous  and  is  contained  in  a  Fabry- 
Perot  cavity. 

3.  For  each  vibrational  level,  the  rotational  levels  have  a  Boltzmann 
distribution  at  the  translational  temperature. 

4.  Only  one  transition  within  a  band  can  lase  at  a  given  time.  This 
is  the  transition  with  maximum  gain  in  the  band  and  is  always  in 
the  P-branch.  During  lasing  on  this  line,  the  gain  at  line  center 
is  held  constant  at  its  threshold  value. 

5.  Initiation  is  modeled  by  the  instantaneous  introduction  of  a  finite 
concentration  of  F  atoms  into  an  initial  gas  mixture  of  F2>  H2> 
and  Ar  diluent,  as  might  be  caused  by  a  very  fast  flashlamp  or 
electrical  discharge.  A  finite  time  history  of  F2  dissociation 
such  as  that  produced  by  most  initiation  methods  can  be  modeled 
and,  indeed,  might  be  of  interest  for  modifying  laser  performance. 
Because  of  the  additional  parameters  involved,  such  modeling  was 
deemed  outside  the  scope  of  the  present  study. 
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With  similar  assumptions  and  a  simpler  form  of  the  reaction  system, 
an  analytical  solution  for  a  two-level  model  of  a  continuous  chemical  laser 
was  obtained.  ^  The  present  study  extends  that  model  to  its  multi-level 
equivalent  and  considers  numerical  solutions  for  a  pulsed  HF  laser.  The 
simulation  is  used  to  investigate  the  effect  of  the  optical  and  chemical  parame¬ 
ters  on  the  laser's  performance. 
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II.  THE  CHEMICAL  LASER  COMPUTER  SIMULATION 


This  section  describes  briefly  the  formulation  of  the  chemical  laser 
computer  simulation.  A  detailed  description  is  given  in  Ref.  9. 

The  chemical  reactions  may  be  written  as 

kfr 

Ea  .  N.  Z==zY'  P  .  N.  (i) 

n  1  ,  Lmt  ri  i 

i  kbr  i 

where  N.  is  the  molar  concentration  of  species  i,  a  .  and  p  .  are  stoichio- 
1  r  ri  rn 

metric  coefficients,  and  and  k^  are  forward  and  backward  rate  coeffi¬ 
cients.  (The  Nomenclature  gives  detailed  definitions  for  symbols  used  in  the 
text.  ) 

The  rate  of  change  of  concentration  N(v)  for  vibrational  level  v  is 
given  by 


*eh<v>  '*■  <v-'’  JL>  <2> 

where  the  +ra{j  terms  are  rates  of  concentration  change  resulting  from  lasing 
into  and  out  of  level  v.  Lower  level  rotational  quantum  numbers  J  and 
are  selected  as  those  giving  maximum  gain  for  transitions  v  +  i  — v  and 
v  —v  -  1,  respectively.  Chemical  reactions  cause  a  rate  of  change  i|» 
which  is 


+ch(vl  ■  E'Prl  -  ■  W 

r 


(3) 


where 


fr 


Nj«rj 


Preceding  page  blank 
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and 


Jbr 


=  k 


br 


IIn.p  . 
i  j  rJ 


The  laser  cavity  is  assumed  to  have  a  uniform  photon  flux  with  active 
medium  length  L  and  mirror  reflectivities  Rq  and  R^.  Lasing  initiates  for 
any  v  +  1  -*v  band  when  the  gain  a  for  the  highest  gain  transition  of  that  band 
is  equal  to  the  threshold  gain  a  where 

®th  =  “  2L  *n<RoRL^  ^ 


The  gain  is  assumed  to  be  constant  over  length  L.  Only  P-branch  transitions 
need  be  considered;  thus,  the  gain  is  given  by  Ref.  9 


o(v,J)  = 


hN  . 
_ / 

4ir 


^c(v,  J)  0  (v,  J)  B  (v,  J) 


r  2j  + 1 

[2J  -  1 


N(v  +  1,  J  -  1)  -  N(v,  J)  (5) 


where  v  and  J  refer  to  the  lower  level  of  the  transition.  The  wave  number 
of  the  transition  is  wc(v,  J),  and  B(v,  J)  is  the  Einstein  isotropic  absorption 
coefficient  based  on  the  intensity.  Line-broadening  constants  and  resonance 
constants  used  in  the  Voigt  profile  at  line  center  0(v,  J)  are  those  of  Ref.  9. 

A  Boltzmann  distribution  at  the  translational  temperature  is  taken  for 
the  rotational  populations,  hence 


N(v,  J)  =  N(v)  <2J  * 

Qr  (T) 


-hcEj/kT 


(6) 


V  V 

where  values  of  rotational  partition  function  Q  (T)  and  rotational  energy  E. 

10  r  J 

are  from  data  of  Mann,  et  al. 
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The  energy  equation  for  a  constant  density  gas  is  written  in  the  form 


yNC  *1.  d£  =  _  p  _y^H 

^  i  pi  dt  dt  L  dt  i 


where  C  ^  is  the  specific  heat  at  constant  pressure  and  IL  the  molar  enthalpy 
of  species  i,  p  is  the  pressure,  and  is  the  output  lasing  power  per  unit 
volume.  This  power  is  given  by 


PL|,»'£hcNA<.c(v,J)+ri<i(v,J) 

v=0 


and  the  pulse  energy  released  in  time  t  is 

EL(t]  =  J  PL(t')dt'  (9) 

0 

Chemistry  and  thermodynamics  determine  the  production  of  excited 
species  until,  at  some  time  tp  and  for  some  J  =  Jq,  the  gain  on  a  given 
vibration-rotation  transition  reaches  or  At  this  time,  the  laser  pulse  begins. 
Then  Eqs.  (2),  (4),  and  (7),  along  with  the  sta  e  equation,  are  solved  for  the 
temperature,  pressure,  concentrations,  and  the  J  for  each  band  that  is  lasing 
as  a  function  of  time.  During  lasing,  the  value  of  J  may  change;  this  process 
will  be  referred  to  as  a  "J-shift."  Lasing  terminates  when  all  gains  drop 
below  or 
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HI.  MECHANISMS  OF  THE  PULSED  HF  LASER 


A.  GENERAL  ENERGY  INTERPRETATION 

The  initial  gas  mixture  in  the  laser  cavity  has  a  given  concentration 
ratio  F^i^iFiAr,  where  the  F  atoms  represent  varying  degrees  of  F^  dis¬ 
sociation  by  some  initiation  scheme.  The  diluent  Ar  is  used  mainly  for  tem¬ 
perature  control.  A  typical  time  history  of  the  total  pulse  power  is  shown  in 
Fig.  1,  and  curves  depicting  rate  of  laser  energy  release  for  the  bands  lasing 
are  shown  in  Fig.  2(e).  For  this  case,  the  initial  F^rH^tFiAr  is  1:1:0.1:50, 
where  one  unit  of  concentration  is  equivalent  to  9.4072  X  10*^  mole/cc.  This 
unit  is  used  throughout  this  work.  Initial  temperature  and  pressure  are  300*  K 
and  1.207  atm,  respectively.  Mirror  reflectivities  are  0.8  and  1.0,  and  the 
active  medium  length  is  100  cm.  This  case  serves  as  a  standard,  and  excur¬ 
sions  from  it  comprise  the  parametric  study. 

In  the  course  of  the  reaction,  the  temperature  and  the  concentrations  of 
HF(v)  increase,  causing  the  J  value  for  maximum  gain  on  a  given  vibrational 
band  to  incr^asi..  In  the  calculation  when  lasing  shifts  to  the  next  higher  J, 
there  is  a  discontinuous,  but  small,  increase  in  the  rate  of  photon  emission. 

In  Fig.  2(hj,  the  laser  energy  emission  rate  for  the  2-1  band  is  shown  as  a 
function  of  time  with  all  the  J-shift  details  included.  The  J's  used  to  label 
J-shifts  on  the  2  -1  band  are  associated  with  the  v=l  level.  The  discontinuities 
without  labels  a  rs;  caused  by  J-shifts  on  adjacent  bands.  For  simplicity, 
discontinuities  have  been  smoothed  in  Figs.  1  and  2(a). 

The  following  paragraphs  give  an  interpretation  of  the  major  mechanisms 
of  the  standard  case.  The  initial  emission  "spike"  results  from  the  pumping 
of  HF(v)  levels  by  Reaction  8-10.  ^  During  this  spike,  approximately  80%  of 
the  initial  F-atom  concentration  is  consumed  (Fig.  3).  The  increase  in  the 
rate  of  photon  emission  after  the  initial  spike  is  due  to  an  increase  in 
temperature  and  to  the  increasing  importance  of  Reactions  14-20  as  a 
result  of  the  production  of  H  by  Reactions  7-10  and,  hence,  the  operation 
of  the  chain. 


Preceding  page  blank 
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PULSE  POWER,  x  10 3  W/cm 


TIME,  ftsec 


Figure  1.  Power  Hiatory  of  the  Standard  Case 


LASER  ENERGY  EMISSION  RATE,  x  KT W/cm 


3  0 


Figure  2.  Laser  Power  and  Energy  Emission  Rate 

(a)  Laser  power  from  the  first  six  levels  of  HF.  (b) 
Laser  energy  emission  rate  from  the  2-1  band  for 
the  standard  case.  The  new  value  of  J  for  the  lower 
level  is  shown  for  each  J- shift  location. 
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CONCENTRATION,  x  10"'  mole/cm 


TIME,  ft  sec 


Figure  3.  Time  Histories  of  the  Concentrations  of  H,,  F?,  H, 
and  F  for  the  Standard  Case 
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Termination  of  the  lasing  pulse  is  caused  in  part  by  the  depletion  of 
the  H£  and  (Fig.  3).  In  addition,  as  the  fuel  supply  is  consumed,  the 
pumping  rate  decreases  until  chemical  deactivation  of  HF(v)  forces  the  gain 
below  threshold. 

The  set  of  equations  (2)  may  be  solved  for  emission  on  a  particular 
band  4*ra<j(v,  J)  and  then  put  in  the  form 

E  V  ■+,«!<''• « +  E  +Dv+i.v  <">> 

v'=V+l  v'=v+l 


where  *v#  is  the  rate  of  pumping  into  level  v '  by  the  Reactions  8-10  and  15-20, 
and  D  is  the  chemical  deactivation  rate  of  the  v+l,v  inversion.  In 

general,  the  VV  HF-HF  deactivation  reactions  are  not  very  important  when 
compared  to  VT  deactivation;  hence,  the  major  contribution  to  y  is 

deactivation  of  HF(v+l)  by  VT  reactions  and  VV  reactions  with  H^v'). 

In  fact,  these  are  the  only  contributions  to  D  ,  ,  when  anharmonic  effects 
are  negligible.  Of  course,  Eq.  (10)  is  a  result  of  conservation  of  number 
density  of  HF(v),  and  terms  in  the  equation  may  be  interpreted  as 


[Pumping] 
[Rate  J 

Photon 

Level 

Chemical 

Emission 

+ 

Filling 

+ 

Deactivation 

Rate 

Rate 

Rate 

where  level  filling  is  the  filling  of  vibrational  levels  as  needed  to  maintain  the 
required  gain.  The  relative  sizes  of  the  terms  in  Eq,  (10)  for  the  2-1  band 
are  shown  in  Fig.  4,  where  the  rate  of  photon  emission  on  that  band  is  also 
shown  as  a  solid  line.  Note  that  Eq.  (10)  is  essentially  a  representation  of 
the  distribution  of  pumping  energy  when  anharmonic  effects  are  small. 

Distribution  of  the  pumping  rates  among  the  vibrational  levels  has  a 
significant  effect  on  the  energy  available  for  lasing.  For  example,  the 
production  of  one  HF(6)  molecule  makes  up  to  six  times  more  energy  available 
for  lasing  than  the  production  of  one  HF(1)  molecule.  Pumping  reactions 
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CHEMICAL  AND  PHOTON  PRODUCTION  RATES,  mole  /cm3-sec 


Figure  4.  Relationship  of  Dominant  Mechanisms  for  the  2-1 
Band  of  HF  in  a  Chemical  Laser 

The  solid  curve  represents  the  instantaneous  rate 
of  release  of  lasing  photons  for  the  2-1  band 
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such  as  8  and  9  that  promote  a  total  inversion  are  more  favorable  for  meeting 
threshold  conditions  than  the  more  even  pumping  distribution  of  Reactions 
15-20.  In  regimes  where  threshold  is  easily  attained  and  in  the  absence  of 
deactivation,  changes  that  increase  the  pumping  of  Reactions  15-20  relative  to 
Reactions  8-10,  such  as  H  or  F2  rich  mixtures,  can  be  expected  to  increase 
pulse  energy,  since  they  supply  more  vibrational  energy  per  mole  of  fuel. 

B.  A  DETAILED  INTERPRETATION  OF  MECHANISMS  FOR  THE 

STANDARD  CASE 

This  section  describes  in  detail  the  mechanisms  that  promote  the 
partial  inversion  and  contribute  to  laser  power  from  the  2-1  band.  In  addition, 
we  shall  assess  the  effect  of  the  photon  flux  on  the  chemical  mechanisms 
associated  with  this  band. 

The  concentration  difference  HF(2)-HF(1)  is  nearly  proportional  to  the 
gain  on  the  2-1  line,  Eq.  (5).  The  relative  importance  of  the  chemical  reac¬ 
tions  and  photon  emissions  that  alter  the  v  =  2  to  v  =  1  inversion  is  depicted  in 
Fig.  5(a).  The  rate  of  mass  transport  by  a  given  reaction  is  L^r  -  L,  A 
weighting  factor  i>r  is  used  to  give  the  appropriate  contribution  of  each  reaction 
to  the  increase  of  HF(2)-HF(1).  For  example,  the  VT  reaction 

HF(2)  +  M-HF(l)  +  M 

decreases  HF(2)  and  increases  HF(1).  The  effective  mass  transported  is  twice 
(L,jr  -  L^)  for  the  reaction;  hence,  vr  =  2.  The  major  contribution  to  mass 
transport  during  the  initial  stages  of  the  reaction  results  from  lasing  on  the 
2-1  and  adjacent  bands.  At  t  =  0,  only  the  pumping  Reactions  8  and  9  are 
present.  In  the  figure,  the  first  bar  shown  is  for  t  =  0+.  The  initial  spike 
in  the  photon  emission  rate  is  a  result  of  the  high  initial  F  concentration, 
and  the  remainder  of  the  curve  (Fig.  1)  may  be  interpreted  by  noting  the  effect 
of  temperature  rise  and  action  of  the  chain,  as  will  be  shown  later. 
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16 


Figure  5.  Production  and  Reduction  Mechanisms  for  the  Gain  of  the  2-1  Line 
(a)  Standard  case,  (b)  Standard  case  with  radiation  suppressed. 


In  the  preceding  section,  the  significance  of  pumping  by  Reactions  16-20 
was  noted.  The  format  of  Fig.  5(a)  tends  to  mask  this  significance,  for  at 
first  glance  Reactions  8  and  9  appear  to  be  the  dominant  chemical  pumping 
mechanisms.  However,  Reactions  16-20  pump  the  upper  vibrational  levels 
and  contribute  to  the  pumping  of  v=2  by  cascading.  We  note  that  the  major 
pumping  and  deactivation  of  the  2-1  inversion  arises  from  lasing  on  the  bands 
1-0,  2-1,  and  3-2.  Lasing  terminates  before  the  fuel  supply  (F^  and  H,)  is 
consumed  because  the  VT  chemical  deactivation 

HF(v)  +  M  —  HF(v  -  1)  +  M  M  =  F,  HF 

becomes  large  enough  to  force  the  gain  below  threshold.  The  rate  of  mass 
transport  by  all  other  reactions  is  small. 

We  examine  the  standard  case  with  radiation  suppressed  by  operating 
the  code  with  reflectivities  set  equal  to  zero.  Production  and  reduction  mech¬ 
anisms  of  the  2-1  inversion  for  this  zero-power  case  are  shown  in  Fig.  5(b). 
Comparing  Figs.  5(a)  and  5(b),  we  note  that  in  the  absence  of  a  photon  flux 
both  pumping  of  the  HF(v)  levels  and  their  deactivation  rates  are  larger. 
Pumping  is  larger  because  the  gas  is  hotter  (hence,  the  reactions  are  faster), 
as  is  shown  in  Fig.  6,  where  the  temperature  history  of  the  standard  case 
and  the  zero-power  calculations  are  compared.  The  increase  in  deactivation 
is  caused  by  the  presence  of  higher  concentrations  of  HF(v).  The  higher  popu¬ 
lations  in  the  HF(v)  levels  for  the  zero-power  case  also  result  in  an  increased 
significance  of  VV  reactions  (47,  48,  65,  66,  and  68)  of  HF(v)  with  HF(v')  and 
H^lv").  In  addition,  we  see  in  Fig.  5(b)  that  VT  deactivation  of  HF(v)  by  the 
diluent  (Reactions  39  and  40)  begins  to  appear. 

For  the  lasing  case,  the  gain  reaches  threshold  and  remains  at  that 
value,  but  in  a  zero-power  calculation  the  maximum  gains  for  the  lower 
transitions  reach  values  nearly  3  orders  of  magnitude  larger  than  the  threshold 
value  (Fig.  7).  In  the  figure,  locations  where  the  gain  for  a  given  band  becomes 
larger  for  a  new  value  of  J  are  denoted  by  the  new  J.  A  comparison  of  the 
gain  on  the  4-3  band  in  Fig.  7  with  the  4-3  curve  in  Fig.  2(a)  shows  the 
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TEMPERATURE, 


Figure  6.  Temperature  History  of  the  Standard  Case,  with  and 
without  Radiation  Suppressed 
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importance  of  radiative  cascading  and  the  limited  usefulness  of  zero-power 
gain  calculations  for  a  multi-level  system  such  as  this. 

In  Fig.  8,  production  and  reduction  mechanisms  for  the  F-atom  popu¬ 
lation  are  shown.  After  a  transient  period  of  consumption  of  F  atoms  by 
Reactions  7-10,  the  production  and  reduction  mechanisms  are  nearly  the 
same,  which  illustrates  the  presence  of  a  quasi-steady  chain. 

As  indicated  earlier,  the  second  hump  in  the  laser  photon  emission  curve 
may  be  explained  by  noting  the  effect  of  the  temperature  increase.  Rate 
coefficients  characteristic  of  the  important  chemical  mechanisms  are  depicted 
as  a  function  of  time  for  the  standard  case  in  Fig.  9.  The  pumping  rates 
increase,  while  the  principal  deactivation  rate  decreases  with  time.  Hence, 
temperature  variations  of  rate  coefficients  during  the  pulse  tend  to  increase  the 
instantaneous  rate  of  photon  emission  and  reduce  chemical  deactivation  losses. 
At  still  higher  temperatures,  the  VT  deactivation  rate  of  HF  begins  to  increase 
with  temperature,  as  is  illustrated  for  Reaction  32  in  the  insert  in  Fig.  9. 

The  increase  in  the  photon  emission  curve  (Fig.  2)  between  0.  15  and  1. 0  psec 
may  be  attributed  to  an  increase  in  pumping  with  the  rise  in  temperature. 

This  may  be  seen  by  comparing  Fig.  9  with  Fig.  3.  During  this  time  interval, 
the  reactant  concentrations  decrease  slightly,  whereas,  due  to  the  temperature 
rise,  the  pumping  reaction  rates  increase  by  a  factor  of  2. 

From  the  previous  discussion,  it  is  clear  that  the  important  rate 
coefficients  give  a  more  favorable  competition  between  stimulated  emission 
and  deactivation  with  a  temperature  somewhat  above  300*  K.  With  pumping 
increased  relative  to  deactivation,  the  pulse  energy  is  larger  because  more  of 
the  fuel  is  consumed  while  the  reacting  mixture  is  above  threshold.  The 
effect  of  a  temperature  increase  on  the  pulse  duration  is  less  certain.  While 
the  overall  reaction  is  faster,  the  increased  pumping  relative  to  deactivation 
tends  to  hold  the  mixture  above  threshold  for  a  longer  time. 
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Figure  8.  Mechanisms  Producing  and  Reducing  F-atom 
Concentration  for  the  Standard  Case 
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FORWARO  RATE  COEFFICIENTS,  x  I0'2  cmVmole-sec 


(32)  HF(2)  +  M  —  HF(1)  +  M 
(9)  F  +  H2  (0)  -  HF (2)  +  H 

(16)  H  +  F2  -  HF(2)+  F 


*32  I012 


TIME,  /xsec 


Figure  9.  Time  Histories  of  Important  Rate  Coefficients  for  the 
Standard  Case 
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IV.  PARAMETRIC  VARIATION 


The  laser  pulse  may  be  characterized  in  part  by  the  pulse  duration  t^ 
and  total  pulse  energy  E^(tj).  In  the  parametric  study,  we  encounter  many 
pulses  with  long  low-power  tails;  hence,  the  time  t required  to  release 
one-half  the  pulse  energy  E^(t^)/2  is  also  a  useful  characteristic.  The  effect 
of  varying  the  initial  chemical  concentrations,  cavity  conditions,  and  rate 
coefficients  is  discussed  by  interpreting  the  variations  in  these  pulse 
characteristics . 

A.  THRESHOLD 

Concentrations  and  temperature  are  given  the  same  initial  values  as  in 
the  standard  case,  while  the  characteristics  of  the  optical  cavity  (mirror 
reflectivities  and  cavity  length)  are  varied.  This  is  equivalent  to  varying  c*th 
as  can  be  seen  from  Eq.  (4).  In  Fig.  10,  the  variation  of  pulse  energy  with  ath 
is  shown.  Additional  abscissa  scales  for  equivalent  values  of  and  L  are 
also  shown.  In  Fig.  10,  the  pulse  energy  is  zero  for  >1.0.  This  is 
consistent  with  expectations  based  on  the  maximum  gain  shown  in  Fig.  7. 

As  a^  approaches  zero,  the  pulse  energy  approaches  a  maximum;  of  course, 
practical  limitations  on  cavity  losses  and  cavity  length  can  make  realization 
of  this  maximum  difficult.  One  aspect  of  decreasing  a^  is  that,  in  addition 
to  increasing  E^(tj),  the  pulse  length  characterized  by  tj^  or  tj  i®  also 
extended,  as  shown  in  Fig.  11.  For  decreasing  values  of  a^,  the  pronounced 
increase  in  tj  relative  to  tj^  *s  resu^  a  l°ng  low-power  tail.  The  time 
required  for  release  of  the  majority  of  the  pulse  energy  is  relatively  insensitive 
to  a^«  as  seen  by  the  small  variation  of  t 

Since  more  energy  is  extracted  for  low  a^>  the  temperature  of  the 
reacting  mixture  is  also  lower.  As  a  ^  is  decreased,  the  vibrational  popula¬ 
tions  are  smaller  and  the  rate  of  VT  chemical  deactivation  is  reduced,  even 
though  the  rate  coefficients  for  these  processes  are  higher  because  of  the 
lower  temperature.  Examination  of  the  reaction  rates  for  the  2-1  band  shows 
that  for  smaller  values  of  a^  the  rate  of  pumping  is  larger,  the  instantaneous 
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Figure  10.  Effect  of  Threshold  Gain  on  Pulse  Energy 
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Rq  (Rl=  1.0,  L  =  100 cm) 
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Figure  11.  Effect  of  Threshold  Gain  on  and 
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power  is  larger,  and  the  rate  of  level  filling  (energy  storage)  becomes  smaller. 
The  mixture  reaches  threshold  at  an  earlier  time  and,  thus,  the  F-atom  con¬ 
centration  is  larger,  and  pumping  at  pulse  initiation  is  higher  and  continues  to 
be  so  throughout  the  run. 

B.  DILUENTS 

Diluents  contribute  to  the  heat  capacity  of  the  total  system  and  act  as 

catalysts  for  the  VT  deactivation  reactions  and  recombination  reactions. 

5 

Since  the  VT  deactivation  rate  of  HF(v)  with  Ar  is  about  10  times  smaller 
at  300 °K  than  the  VT  rate  of  HF(v)  with  itself,  Ar  contributes  little  to  chem¬ 
ical  deactivation;  its  major  role  is  temperature  control. 

The  effect  of  variations  in  dilution  is  studied  for  two  initial  concentra¬ 
tions  of  F.  For  initial  F  equal  to  10%  of  initial  F^,  the  following  trends  are 
noted  as  the  concentration  of  Ar  is  increased. 

1.  The  temperature  rise  throughout  the  reaction  is  less, 
as  expected. 

2.  The  laser  power  decreases  and  the  second  peak  in  the  power 
curve  disappears.  Also,  pumping  and  level  filling  rates 
decrease  while  deactivation  rates  increase  as  a  result  of  the 
lower  temperature. 

3.  The  pulse  energy  decreases  gradually  (Fig.  12)  because  of 
the  less  favorable  competition  described  in  Item  2. 

4.  The  pulse  duration  tj  and  time  parameter  t^^  increasei 
as  expected,  for  lower  temperatures  (Fig.  13). 

Initial  F  concentrations  in  this  10%  range  and  above  emphasize  Reactions 
8-10  relative  to  Reactions  15-20. 

For  initial  F  equal  to  1%  of  initial  F^,  pumping  rates  throughout  the 
reaction  are  much  lower.  Deactivation  is  more  significant  and  the  pulse 
energy  is  lower.  As  the  concentration  of  Ar  is  increased,  the  temperature 
throughout  the  pulse  is  lower  and  the  pulse  length  is  increased.  Further 
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Figure  12.  Effect  of  Diluent  Concentration  on  Pulse  Energy  with 
Two  Different  Initial  F-atom  Concentrations 
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(a)  Effect  on  tj/2  with  two  initial  F-atom  concentra 
tions.  (b)  Effect  on  pulse  length  tj  with  two  initial 
F-atom  concentrations. 


increases  in  Ar  lower  the  temperature  sufficiently  such  that  the  stimulated 
emission-deactivation  competition  forces  the  pulse  length  to  decrease  and  a 
maximum  pulse  duration  is  found. 

C.  INITIAL  F?  DISSOCIATION 

The  initial  F-atom  concentration  is  increased  while  the  stoichiometric 
balance  is  maintained  between  hydrogen  and  fluorine.  Other  initial  conditions 
are  held  constant.  As  the  concentration  of  F  is  increased  from  very  low 
values,  the  following  trends  are  noted. 

1.  The  rate  of  pumping  increases,  laser  power  increases,  and  the 
rate  of  level  filling  increases. 


2. 

3. 


The  faster  pumping  rates  result  in  smaller  pulse  durations  t^ 
and  smaller  values  of  (Fig.  14). 


The  pulse  energy  increases  until  the  mass  fraction  of  F 
reaches  0.2  (Fig.  15).  Further  increases  in  F  result  in  a 
dominance  of  reactions  7-10  throughout  the  pulse.  The  second 


peak  in  the  power  curve  disappears,  and  E^(tj)  and  Ev+^  v  for 
each  band  decrease.  Finally,  for  total  dissociation  of  F^, 
Reactions  14-20  do  not  contribute  to  the  pulse  energy  and  E^j, 


E54*  and  E65 


are  zero. 


Since  the  pumping  Reactions  7-10  release  only  about  one-third 
as  much  energy  as  Reactions  14-20,  the  reaction  temperature  is 
lower  for  large  initial  F.  Such  behavior  leads  to  the  stimulated 
emission-deactivation  competition  suggested  in  the  discussion 
of  the  standard  case . 


5.  For  small  concentrations  of  F,  there  is  a  strong  competition 
between  stimulated  emission  and  deactivation.  As  initial  F^ 
dissociation  increases,  the  relative  importance  of  VT  deactiva¬ 
tion  decreases  until  Reactions  23-30  make  a  significant  contribution. 
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Figure  14.  Effect  of  Initial  Dissociation  on  Pulse  Length 
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Figure  15.  Effect  of  Initial  Dissociation  on  Pulse  Energy 
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For  the  larger  Initial  F  concentrations,  the  mixture  is  below 
stoichiometric  for  H^.  To  assess  this  effect,  the  same  variations  of  F  were 
made  with  the  initial  concentration  doubled.  The  results  of  these  compu¬ 
tations  are  shown  in  Fig.  15  as  a  dashed  curve.  There  is  a  marked  increase 
in  pulse  energy  for  high  initial  F.  There  is  also  a  slight  decrease  in  tj  and 
tj/2  as  a  result  of  the  increased  reaction  rates. 

D.  INITIAL  TEMPERATURE 

The  initial  temperature  was  varied  while  all  other  conditions  were 
maintained  at  the  values  of  the  standard  case.  As  T^  increases,  the  tempera¬ 
ture  throughout  the  pulse  is  higher,  reaction  rates  are  faster,  and  characteris¬ 
tic  times  become  shorter  (Fig.  16).  For  T.  above  400*K,  the  hump  in  the 
power  curve  disappears;  i.e.,  throughout  the  pulse  the  rate  of  pumping 
decreases  monotonically,  since  the  effect  of  temperature  increase  is  less 
significant.  In  addition,  the  pulse  energy  begins  to  decrease  slightly;  how¬ 
ever,  the  overall  effect  of  T.  on  pulse  energy  is  small. 

E.  INITIAL  H-.  CONCENTRATION 

The  initial  concentration  of  is  varied  while  other  initial  conditions 
are  those  pf  the  standard  case.  The  pulse  energy  and  pulse  duration  are  both 
quite  sensitive  to  the  initial  concentration  for  I^-lean  mixtures  (Fig.  17). 

In  this  regime,  increases  in  are  accompanied  by  increased  pumping, 
increased  pulse  energy,  and  increased  pulse  duration. 

The  effect  of  variations  in  is  very  small  for  H^-rich  mixtures.  As 
the  concentration  is  increased,  the  pumping  rate  increases,  the  pulse 
energy  increases,  and  the  pulse  duration  decreases  because  of  the  faster 
reaction  rate.  Note  that  the  time  tj^  required  to  deliver  one-half  the  pulse 
energy  is  virtually  independent  of  initial  H^.  Of  course,  very  high  con¬ 
centrations  give  rise  to  VV  and  VT  deactivation  of  HF(v)  with  H^fv'  ). 

F.  VARIATION  OF  IMPORTANT  RATE  CONSTANTS 

The  effect  of  variations  in  dominant  rate  constants  was  studied  with 
variations  from  the  standard  (see  Appendix)  selected  to  include  the  major 
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uncertainties  in  the  rates  at  the  time  of  the  study.  The  effect  of  these  changes 

on  the  pulse  characteristics  and  the  energy  released  from  the  band  v  +  1  -*  v, 

E  . ,  may  be  seen  in  Table  I.  The  variations  in  the  rates  studied  are 
v+l,v,  7 

described  below. 


Rate  1: 


Rate  2: 


Rate  3: 


Rate  4: 


Rate  5: 


Pumping  into  vibrational  levels  i-6  is  reduced  by  a  factor  of  2, 
while  the  total  reaction  rate  is  held  constant  by  increasing  the 
rate  of  production  of  HF(0).  As  expected,  the  pulse  energy  and 
the  pulse  length  decreased  with  the  lower  pumping.  Also,  the 
energy  in  the  1-0  band  is  much  less  because  the  rate  of  popula¬ 
ting  HF(0)  is  much  higher. 

The  rates  for  HF-HF  VV  reactions  are  increased  by  a  factor 
of  100.  This  change  has  a  very  small  effect  on  the  pulse  energy 
and  the  characteristic  times;  however,  it  does  redistribute  the 
energy  emitted  by  the  various  bands. 

The  rates  for  HF-H^  VV  reactions  are  increased  by  a  factor  of 
100.  In  this  case,  the  pulse  energy  is  essentially  constant,  the 
energy  in  the  2-1  band  is  higher  while  the  energy  in  the  1-0  band 
is  lower.  Reaction  65  becomes  a  significant  deactivation  mech¬ 
anism  for  the  1-0  inversion,  and  for  this  mixture  the  combination 
of  Reactions  65  and  66  aids  the  2-1  inversion.  In  addition,  since 
the  1-0  band  contributes  energy  early  in  the  pulse,  t  j  ^  increases 
slightly. 

The  changes  in  Rates  1,  2,  and  3  are  combined.  A  much  lower 
pulse  energy  and  shorter  pulse  duration  are  observed,  as  expected. 

The  backward  rates  for  F  +  H2^HF(v)+  H,  v=  4,  5,  6  are  reduced 
by  a  factor  of  100.  The  pulse  energy  increases  10%.  A  marked  in¬ 
crease  in  the  energy  released  in  the  bands  6-5,  5-4,  and  4-3  is 
noted  as  a  result  of  the  rate  change,  and  the  energy  released  in 
the  lower  bands  is  higher  as  a  result  of  cascading. 
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Table  I.  Effect  of  Rate  Changes  on  Pulse  Characteristics 
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V.  CHEMICAL  EFFICIENCY 


Of  the  31.56  kcal/mole  released  by  Reactions  7-10,  22.6  kcal/mole  is 
supplied  as  vibrational  energy  of  HF(v)  with  the  pumping  rates  listed  in  the 
Appendix.  Reactions  14-20  release  98.04  kcal/mole  of  which  43.1  kcal  of 
vibrational  energy  is  available.  By  combining  these  reactions,  we  find  50.7% 
of  the  chemical  energy  converted  to  vibrational  energy.  After  the  threshold 
condition  and  chemical  deactivation  are  considered,  only  a  fraction  of  this 
vibrational  energy  is  available  for  laser  emission. 

A  thermodynamic  efficiency  of  the  laser  system  would  account  for  the 
energy  required  to  dissociate  F^  upon  pulse  initiation.  However,  for  simplic¬ 
ity,  efficiency  q  will  be  taker  as 

E^(tj)  X  100/mole  F2 
11  '  129.6  kcal/mole  F^ 

Thus,  pumping  rates  alone  require  q  to  be  less  than  50.7%  for  the  chain. 

For  the  standard  case,  q  is  20.6%.  Optimization  of  q  has  not  been  attempted. 
However,  we  note  that  the  effects  of  the  parameters  studied  in  this  report 
indicate  the  probable  direction  to  be  taken  for  optimization,  and  that  the 
variation  of  q  is  approximately  like  that  of  EjJtj). 
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VI.  SUMMARY  AND  CONCLUSIONS 


A  computer  simulation  is  used  to  study  the  interaction  of  the  chemistry 
and  the  optical  cavity  of  a  pulsed  HF  laser.  A  detailed  study  of  the  effect  of 
the  optical  cavity,  the  initial  chemical  composition,  the  initial  temperature, 
and  variations  in  uncertain  chemical  rate  coefficients  on  the  pulse  length  and 
pulse  energy  is  made.  Variations  about  a  standard  case  are  used  to  study  the 
effect  of  these  parameters. 

Pumping,  vibrational-translational  deactivation  reactions,  temperature 
control,  and  threshold  gain  are  found  to  be  important  parameters  governing 
the  behavior  of  the  laser.  A  competition  between  stimulated  emission  and 
deactivation  is  observed,  and  conditions  where  the  effect  of  deactivation  is 
minimized  are  indicated.  A  decrease  in  the  threshold  gain  is  accompanied 
by  an  increase  in  pulse  energy  and  pulse  duration.  The  photon  flux  is  found 
to  have  a  large  effect  on  the  populating  mechanisms  of  the  vibrational  levels. 

Argon  is  used  as  a  diluent  for  temperature  control.  An  optimum  con¬ 
centration  of  Ar  is  found  that  maximizes  pulse  energy,  and  at  low  concentra¬ 
tions  of  initial  F,  an  optimum  concentration  of  Ar  is  found  for  maximizing 
pulse  duration. 

Pulse  duration  decreases  with  increasing  initial  dissociation,  while 
an  optimum  dissociation  of  F^  is  found  that  maximizes  the  pulse  energy. 

The  effect  of  varying  the  initial  temperature  of  the  reactants  on  the  pulse 
energy  is  small.  An  increase  in  the  initial  temperature  decreases  the  pulse 
duration,  as  expected. 

An  H£-rich  mixture  is  desirable  to  provide  an  increase  in  pulse  energy 
and  pulse  duration.  The  effect  of  varying  the  initial  H ^  concentration  for 
I^-rich  mixtures  is  small. 

The  variations  of  the  more  uncertain  rate  coefficients  at  the  time  of 
this  study  reveal  that  the  pulse  characteristics  are  sensitive  to  changes  in 
the  total  pumping  reaction  rates  and  to  the  distribution  of  pumping  in  the 
vibrational  levels.  These  characteristics  are  not  very  sensitive  to  changes  in 
the  VV  rates  of  HF(v)  with  HF(v' )  or  ^(v' ' ).  The  major  effect  of  these  change 
is  redistribution  of  energy  emitted  by  the  lasing  bands. 

Preceding  page  blank 
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APPENDIX 


THE  REACTION  SYSTEM 


The  chemical  kinetic  model  and  corresponding  rate  coefficients  used 

in  this  study  are  listed  in  the  following  table.  Chemical  rate  coefficients  for 

the  HF  laser  system  are  from  a  compilation  prepared  by  N.  Cohen  of  our 

laboratory  in  December  1970.  Subsequently,  Cohen  has  reported  a  more  final 

12 

review  of  these  rates.  Implications  of  the  difference  in  rates  between  the 
table  and  Ref.  12  are  not  serious  and  are  covered  in  the  discussion.  Some 
of  the  catalytic  species  listed  at  the  bottom  of  the  table  are  multiplied  by  a 
constant.  This  is  equivalent  to  using  a  rate  coefficient  whose  value  is  larger 
by  this  factor.  The  thermodynamic  data  are  taken  from  the  JANAF  (Joint- 
Army-Navy- Air  Force)  publication  distributed  by  Dow  Chemical  Corporation. 
However,  the  thermodynamic  data  of  excited  species  are  generated  at  this 
laboratory. 


Table  of  Reactions  and  Rate  Coefficients 


^Catalytic  apecias: 

M  *  H,  F,  Ar,  HF(0) . HF(8),  4*H2(0),  4*H.U),  4*H2<2),  F2 

M,  =  20*H,  F,  Ar.  HF<0) . HF<8|,  2.  2. 2.5*H,(2>,  F, 

M4  1  F  c  c 

Mj  -  H,  F,  Ar,  HF(0) . HF{8),  H,(0),  H,U>.  H,(2?,  F, 

M5  =  H,  Ar,  H.IO*.  H»<2),  F, 

Mi  *  HF(O),  ....  HF(8> 

> 

Th«  rat*  coefficient!  k+  and  k_  dcaignat*  forward  and  backward  rataa,  respectively , 
with  unita  in  tarma  of  molea,  cm\  and  acc.  For  aach  raactlon,  th*  miaaing  rata 
coefficient  ia  determined  from  tha  equilibrium  conatant. 

^Tha  quantity  0  *  -(toJ/RT),  whar*  th*  temparatur*  T  ia  in  *K  and  R  ia  1.987  cai/mol*-*K. 
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